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^3 ■ ABSTRACT 

9-i The growth of supermassive black holes appears to be driven by both galaxy mergers 

and 'secular' processes that occur in their absence. In order to quantify the effects of 
secular evolution on black hole growth, we require a sample of active galactic nuclei 
(AGN) in galaxies that have formed without significant mergers, a population that 
heretofore has been difficult to locate. Here we present an initial sample of 13 AGN 
in massive (M* > 10 10 M Q ) bulgeless galaxies — which lack the classical bulges 
believed inevitably to result from mergers — selected from the Sloan Digital Sky 
Survey using visual classifications from Galaxy Zoo. Parametric morphological fitting 
confirms the host galaxies lack classical bulges; any contributions from pscudobulgcs 
are very small (typically < 5%). This is the largest such sample yet assembled. We 
compute black hole masses for the two broad-line objects in the sample (4.2 x 10 6 
and 1.2 x 10 7 M©) and place lower limits on black hole masses for the remaining 
sample (typically A/bh It 10 6 M©), showing that significant black hole growth must 
be possible in the absence of mergers. 

The black hole masses are systematically higher than expected from established 
bulge-black hole relations. However, if the mean Eddington ratio of the systems with 
measured black hole masses (L/L-Edd ~ 0.065) is typical, 10 of 13 sources are consistent 
with the correlation between black hole mass and total stellar mass. That pure disk 
galaxies and their central black holes may be consistent with a relation derived from 
elliptical and bulge-dominated galaxies with very different formation histories implies 
the details of stellar galaxy evolution and dynamics may not be fundamental to the 
co-evolution of galaxies and black holes. 
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1 INTRODUCTION 



* This publication has been made possible by the participation 
of more than 200,000 volunteers in the Galaxy Zoo project, 
f E-mail: brookc.simmons@yale.edu 
X Einstein Fellow 
§ www.sepnet.ac.uk 



Constraining the contribution of mergers to the evolution of 
the galaxy population is one of the fundamental challenges 
in modern galaxy formation theory. Galaxies have long 
been believed to form hierarchically, building up to their 
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observed sizes through a ser ies of mergers (|White fc Reesl 
1 19781 ; iKauffmann et all 1 19931 ). The merger history of each 
galaxy thus contributes significantly to the galaxy's stel- 
lar and gas dynamics, and is also thought to drive the co- 
cvolution of a gala xy with its cen t ral supermassive black 
hole (SMBH; ISanders et all 1 19881: |Pi Matteo et al.1 120051 : 



ICroton et aTlbOOrj ; IHopkins et al.H2006l , 120081 ). Given such 
a wide ranging effect, distinguishing between the effect of 
mergers and that of dynamically cold processes such as stel- 
lar feedback and cold accretion (collectively known as 'sec- 
ular' evolution) is difficult. In this paper we present a sam- 
ple of massive galaxies which cannot have had a significant 
merger in their history, discuss their properties, and demon- 
strate for the first time that significant black hole growth 
(to Mbh ^ 10 6 M©) is possible without the advent of a 
significant merger. 

A galaxy's morphology contains signatures of its merger 
history. In particular, the assembly of massive disk galaxies 
through m ergers inevitably produces a central bulge compo- 
nent fe.g..lToomrdll977l;IWalker et al.lll99tj ; IHopkins etal] 
l201ll ; iMartig et all l2012h . The bulge is dynamically hot, 
rising vertically above the disk, and has a steeper den sity 
profile than an exponential disk (|de Vaucouleurdll953h . A 
galaxy lacking a central bulge thus has a formation history 
free of significant mergers, with a strong limit on the mass 
ratio between the main galaxy and any accreting satellite 



galaxie s (~ 1 : 10; IWalker et al.|[l996l ; IHopkins et ai1l201ll . 



though I Brook et al.l 20121 suggest the ratio may be slightly 
higher) . 

Such bulgeless galaxies, with a purely secular formation 
history, might be expected to be rare in a hierarchical sce- 
nario. The presence amongst the galaxy population of large 
bulgeles s galaxies thus present s a serious challenge to this 
picture (jKormendv et al]|2010i ). as they cannot have under- 
gone a significant merger yet have assembled stellar masses 
of M* > 10 10 M Q . 

Additionally, the well-establishe d co-evolution of galax- 
ies and their central SMBHs (e.g., iMagorrian et al.l ll99Sl : 
iKormendv fe Gebhardtl l200ll ) shows a strong correlation 
between a host galaxy's bulge properties and the mass 
of its central black hole (e.g.. iFerrarese fc Merrittl l2000l; 
Tremaine et al.ll2002i ; [Marconi fc Huntll2003l ; lHaring fc Rrxl 
This has led to the prevalence of major- 



20041 ). 



merger-driven theori e s for black hole - galax y co-evolution 
JSanders et al.l ll98S) ; |Pi Matteo et al.1 120051 ; ICroton et all 
120061 ; IHopkins et all I200SI ), but a growing body of recent 
work suggests minor mergers and secular processes may be 
a mo re typical means o f growing a SMB H, both locally 



growing black holes. Selecting galaxies that lack classical 
bulges (as opposed to galaxies with a more varied history 
of both secular and merger-driven evolution) enables the 
isolated study of black hole growth in the absence of merg- 
ers. These galaxies provide a strong challenge to models of 
galaxy formation, requiring substantial and ongoing secular 
growth of a central black hole. 

We aim to use this rare population to assess whether 
these galaxies fall on the same galaxy-SMBH relations seen 
in galaxies with more merger-driven histories. By compar- 
ing upper limits on bulge masses to black hole masses from 
broad emission lines or lower limits on black hole masses 
using Eddington limits, we assess the sizes to which black 
holes can grow over their lifetimes due to secular processes 
alone. 

Section[2]describes the methods used to select bulgeless 
galaxies with growing black holes. Section [3] presents the 
sample of host galaxies, with Section U detailing how the 
black hole masses and lower limits are calculated. In Section 
[5] we discuss how bulgeless AGN host galaxies inform our 
understanding of the co-evolution of black holes and galax- 
ies. Throughout this paper, we assume Ho = 71 km/s/Mpc, 
fijvi = 0.27 and $1a = 0.73, consistent wi th the most recent 
WMAP cosmology (|Komatsu et al.ll201ll ). 



2 FINDING BULGELESS AGN HOST 
GALAXIES IN GALAXY ZOO 

We use visual morpholog ies drawn from the G alaxy Zoo 2 
project, first described in lMasters et al.l |201ll ). to assemble 
a sample of approximately 10, 500 disk galaxies drawn from 
SDSS which appear bulgeless or nearly bulgeless. We then 
select a much smaller sample of bulgeless galaxies which host 
growing supermassive black holes. This section describes 
first the initial selection of disk galaxies, followed by the 
AGN identification and hence a conservative selection of bul- 
geless AGN host galaxies. 

Galaxy Zoo volunteers are asked to classify randomly 
chosen colour images of SDSS systems by clicking buttons 
in response to a set of descriptive questions arranged into a 
decision tree. The most relevant here is a question which asks 
volunteers to classify the bulges of systems already identi- 
fied as face-on spirals into one of four categories : no-bulge, 
just-noticeable, OBVIOUS and dominant. A full descrip- 
tion o f the Galaxy Zoo decision tree is given in lMasters et al.l 
(|201ll ). 



(e.g., I Greene et al. 2010; Jian g et al 
redshift (e.g.. | Simmons et al.l 20 111: 



201 ll) and at higher 



ISchawinski et al.l l201ll . |2012| ; iKocevski et al.l |2012| ). How 

ever, owing in part to the compounded rarity of both mas- 
sive, bulgeless galaxies and active galactic nuclei (AGN), 
the the extent to which a SMBH can grow in the absence of 
merger processes remains difficult to characterise. 

This paper uses morphological classificati ons from the 
Galaxy Zoc0 project l|Lintott et al.ll200Sl . I2O1TI ) to construct 
a sample of bulgeless galaxies drawn fr om the Sloan Digi- 
tal Sky Survey (SPSS; lYork et al.ll2000h that host actively 



Cistcrnas et al.l iS 2>1 AGN BiaS: Point Sources Can Mimic Bul S es 



www.galaxyzoo.org 



Galaxy Zoo obtains many independent classifications of each 
system because it is a citizen science project with a large 
number of participants. While this approach has many ad- 
vantages over classification by either an individual or a small 
group of experts, it is still prone to the biases inherent 
in morphological classification. These are particularly acute 
when dealing with small bulges in systems with AGN, where 
the presence of a nuclear point source can be confused with 
a central bulge. 

In order to investigate the size of this effect on mor- 
phological classification, simulated AGN were added to a 
subsample of images in the most recent iteration of Galaxy 
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Zoo, which uses data drawn from large Hubble Space Tele- 
scope surveys including GOODS, GEMS and COSMOS. 
These AGN host galaxy simulations, creat ed using a similar 
metho d to those described in Section 3.3 of lSimmons fc Urrvl 
l|2008l ). are at higher redshifts than the AGN hosts in the 
SDSS sample discussed here. However, the much higher res- 
olution of the HST Advanced Camera for Surveys (ACS) 
means the simulations cover a similar spatial resolution as 
images of SDSS galaxies in Galaxy Zoo 2: an ACS image of 
a galaxy at z — 1 has the same resolution, in kpc per pixel, 
as a SDSS image of a galaxy at z = 0.053. The results of this 
test of AGN bias on morphological clas sification are thu s di- 
rectly applicable to the present study (|Koss et al.ll201ll also 
find a similar parallel between host galaxy simulations using 
HST images of galaxies at z ~ 1 and lower-redshift SDSS 
images) . 

The synthetic AGN host morphologies show that the 
presence of even a faint nuclear point source can signifi- 
cantly affect the visual classification of a galaxy that would 
otherwise be classified as bulgeless. Among inactive galax- 
ies with a NO-BULGE classification of at least 80%, the ad- 
dition of a nuclear point source with just l/50th the lu- 
minosity of the host galaxy decreases the no-bulge clas- 
sification by at least 50%; those classifications are instead 
transferred mainly to the just-noticeable category, such 
that no-bulge + just-noticeable > 70% for a bulgeless 
host galaxy with a faint simulated AGN. When the nuclear 
point-source luminosity is increased to l/10th that of the 
host galaxy, the NO-BULGE classification decreases to 10% 
or less, with a corresponding and significant increase in the 
just-noticeable and OBVIOUS bulge classifications. As the 
AGN luminosity increases with respect to the host galaxy, 
the visually classified bulge fraction increases substantially. 

While it is well established that parametric morpholo- 
gies can overestimate a bulge contribution if an AGN is 
present but not accoun ted for in a parametric analysis (e.g., 
ISimmons fc Urrvll2008l ). these results from Galaxy Zoo sim- 
ulations demonstrate that this tendency also appears in vi- 
sual classifications of AGN host galaxies. Without care being 
taken to distinguish nuclear activity from a bulge, there- 
fore, a strict selection is likely to reject many truly bulge- 
less galaxies hosting both unobscured and obscured AGN. 
We account for this bias using a combination of a more re- 
laxed initial selection, follow-up visual inspection, and fi- 
nally through parametric separation of host galaxy from 
AGN. 



2.2 Classical Bulges versus Pseudobulges 

One particular difficulty in the discovery and analysis 
of bulgeless galaxies is the distinction between a bulge 
and a pseudobulge. A class ical bulge, as defined by 
iKormendv fc Kennicuttl (|2004l ) . is 'an elliptical living in a 
disk', formed by mergers as discussed above, while a pseu- 
dobulge still retains signs of having been formed from disk- 
driven processes. A pseudobulge is, therefore, unlike a bulge 
in that it is a dynamically cold system, with stars distributed 
in a disk but with a somewhat steeper density profile than 
a typical disk. 

Pseudobulges may also be marked by spiral struc- 
ture or the presence of a nuclear bar or starbursts, and 
may have profiles similar to the exponentials seen in disks. 



Their formation mechanism is a matter o f debate (e.g., 
IKormendv et al.ll2010l ; lFisher fc Drory||2010l ) but, by defini- 
tion, the presence of a pseudobulge in a galaxy is consistent 
with a merger-free history. 

The ability to distinguish between pseudobulges and 
classical bulges is therefore of key importance in determin- 
ing whether a galaxy is truly bulgeless. Because pseudob- 
ulges have light profiles more closely resembling exponential 
disks than classical bulges, parametric morphological fitting 
of a pseudobulge with a Sersic profile should indicate a more 
disk-like profile (where an exponential disk has a Sersic in- 
dex of n = 1) rather than a classical deVaucouleur bulge 
(which has n = 4). A Sersic index of n = 2 is commonly 
used as a divider (|Fisher fc Drory||2008i ) , and this criteria is 
used to distinguish between these categories in Section [2.41 



2.3 Sample Selection 

2.3.1 Broad Parent Sample 

Galaxies were selected where at least 30 people classified 
the system as spiral and face-on, and thus answered the 
question 'How prominent is the central bulge, compared with 
the rest of the galaxy?' We further require that the combined 
classifications no-bulge+just-noticeable ^ 0.7, meaning 
at least 70% of the weighted classifications fall into either of 
these categories. 

Additionally, we re quire that the galaxies are include d 
in the OSSY catalogue (|Oh. Sarzi. Schawinski. fc Yi|[201ll ). 
which provides emission and absor ption line measure ments 
produced by the GANDALF code (|Sarzi et alj|200rj ). This 
produces a parent sample of 10, 488 galaxies with either no 
bulge or a small bulge. Future work will concentrate on re- 
fining this sample to provide an estimate of the population 
density of truly bulgeless galaxies; this paper will concen- 
trate on the subsample which have growing black holes. 



2.3.2 AGN Selection 

Active Galactic Nuclei are selected from the parent sam- 
ple describe d above using the optical line diagn ostic first 
described bv lBaldwin. Phillips, fc Terlevichl (|l98ll ). Requir- 
ing line measurements with signal-to-noise ratio S/N ^ 3 
for each of the [Om] A5007, [Nil] A6583, Ha, and H/3 lines 
produces 5,904 sources from the broad parent sample; their 
positions on the BPT diagram are shown in Figure [T] 

We select galaxi es in the region above both the extreme 
star-formation line of lKewlev ct al. (2001) and the e mpirical 
AGN-LINER separation of ISchawinski et all (|2007l ), which 
are unambiguously AGN hosts; 100 galaxies from the parent 
sample lie in this region of the parameter space. 

Because this AGN selection method requires strong op- 
tical emission from the AGN, all the selected AGN are ex- 
pected to have visually detectable optical point sources. 
As the parent sample was selected using broad bulge- 
classification criteria in order to account for confusion be- 
tween point-sources and small bulges, many AGN hosts se- 
lected will not be truly bulgeless, but will have small bulges. 
To further select a sub-sample of truly bulgeless host galax- 
ies, two authors (BDS and CJL) visually inspected each of 
the AGN+host galaxy images, selecting only those images 
with no indication of an extended bulge regardless of the 
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Figure 1. Emission li ne ratios for u se as a diagnostic of 
AGN activity, followin g I Baldwin et al.l l|l98lh . The solid line 
llKauffmann et al.lliooj) empirically separates pure star forming 
galaxies from composite sources with a com bination of star for - 
mation and AGN activity. The dotted line jKewlev et al.ll20ofl) 
shows the limit for ex treme star formation. The dashed line 
llSchawinski et al.l [20071) shows the empirical AGN- LINER sep- 
aration. Gray points represent galaxies having a summed Galaxy 
Zoo 2 classification of NO-BULGE + JUST-NOTICEABLE ^ 70% and 
emission line S/N ^ 3 for all four emission lines, from which we 
selected 15 sources in the AGN region of the above diagram with 
no visual evidence of a bulge despite an obvious nuclear point 
source. 



point-like nuclear emission. This very conservative cut pro- 
duced a sample of 15 AGN host galaxies that appear to 
be completely lacking a bulge. Figure [5] shows SDSS colour 
images of each system. 

As noted above, this is a lower limit on the AGN 
fraction as optical selection will not find heavily obscured 
sources. We have also excluded galaxies in the transition re- 
gion of the BPT diagram, which will include a subs tantial 
number of AGN (|Trouille. Barger. fc Tremontill201 J ). 

2.4 Parametric Morphological Fitting 



Section 12.11 describes the confusion between point sources 
due to the AGN and any bulge at the centre. In visually 
inspecting our sample we did not exclude sources with a 
point-like source at the centre; for a quantitative decompo- 
sition betwe en galaxy and point so urce parametric fitting 
is necessary (|Simmons fc Urrv||2008h . Parametric fitting has 
the added advantage of enabling both an independent as- 
sessment of host morphology and quantitative separation 
between disk and bulge, providing constraints on possible 
bulge contribution. 

Separation of host galaxy from nuclear emission requires 
careful characterisation of the image point spread function 
(PSF). For ground-based observations such as these, the 
PSF can vary significantly depending on atmospheric con- 



ditions, and is difficult to model purely analytically. We 
therefore use the PSF-creation routines in the IRAF pack- 
age daophot to create a semi-analytical PSF for each image 
based on stars observed near each system. The number of 
stars available for each source varies with the source's dis- 
tance from the galactic plane; between 2 and 40 stars were 
used for each source (the median number of stars used in 
PSF creation is 14). Modelling several stars minimizes noise 
compared to using a single star, whilst still accounting for 
the unique conditions at each epoch of observation. 

We use the two-dimension al parametric fitting program 
GALFIT (|Peng et alj2002l . l201Ch to simultaneously model the 
unresolved nucleus and extended galaxy for each of the 15 
AGN+host galaxies selected above, choosing the r-band im- 
ages for their depth. Initially, we fit each source's central 
regio n with a combination of a single Sersic profile (|Sersid 
Il968h and a central point source. Initial parameters (mag- 
nitude, radius, axis ratio and position angle) were either 
drawn from the catalogue or estimated where they are not 
given by the SDSS catalogue. Initially, the host Sersic index 
is set to n = 2.5 and allowed to vary. This value was chosen 
so as to avoid favouring either an exponential disk (n = 1) 
or a deVaucouleur bulge (n = 4). We find, however, that 
the final best-fit parameters (i.e., with a minimum x 2 ) are 
insensitive to initial indices between 1 ^ n ^ 4, so long as 
the other initial parameters are reasonable. 

The primary purpose of this initial fit is to converge on 
the centroid positions of each component; successive itera- 
tions include the extended regions of the galaxy. In order 
to ensure the extended galaxy is properly fit, we fix the sky 
background to an independently-determined value for each 
individual source. Where present, we also fit and subtract 
nearby bright stars and extended companion galaxies, and 
mask fainter compact sources from the fit. 

Throughout, the primary goal of the lit was to neither 
over- nor under-subtract the galaxy's central region. In most 
cases, this can only be achieved by either masking out asym- 
metric bright features (such as star-forming knots, bars and 
spiral arms) or adding them as additional components of the 
host galaxy fit. We fit these additional features only when 
they are necessary to ensure the disk component is properly 
modelled in the central region of the galaxy. 

Once the single Sersic plus nuclear point-source fit has 
converged to its best-fit solution, we add a small second 
Sersic component with a variable n and initially equally 
bright as the original source, both to constrain the contribu- 
tion of a small extended bulge that may have been visually 
obfuscated by the nuclear point source and to distinguish 
between bulge and pseudobulge. As outlined in Section \2. 21 
compact host galaxy components having light profiles with 
Sersic indices n < 2 are considered pseudobulges, whereas 
components with n > 2 are considered classical bulges. 

Figure [3] shows the SDSS r-band images and the 
residual images after subtraction of the best parametric 
fits. All of the recovered host parameters are reliable be- 
cause the AGN are sig nificantly fainter than their hosts 
l|Simmons fc Urrvl 12008 ). In practice, fitting multiple com- 
ponents to a source increases the uncertainties in recov- 
ered fit parameters over the computa tional uncertainties 
reported by GALFIT (|Peng et alj|20ich . We therefore add 
additional uncertainties to those reported by GALFIT using 
results from extensive parametric host galaxy fitting simula- 
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Figure 2. SDSS colour cutouts of 15 AGN (shown in Figurc[T]as green open circles) with potentially bulgeless host galaxies from visual 
selection. Reading from top left, the images are sorted by ascending redshift, from z = 0.014 to z = 0.19, matching the order in Tablc[T] 
Each cutout is marked at the top left with a scale bar representing 5". 



3 SAMPLE PROPERTIES 

3.1 Are these AGN host galaxies really bulgeless? 

The remaining sample of 13 host galaxies are all well fit by 
a model consisting of a dominant disk and a nuclear point 
source, providing strong constraints on the maximum con- 
tribution of a small bulge component. In 3 cases, we do not 
detect a second, compact host component. For the other sys- 
tems, examination of residuals from fitting only the extended 
disk + nuclear point-source components shows clear signs of 
a small extended component in the center; in all but one 
case this additional component has a Sersic index consistent 
with a pseudobulge (n < 2 within the la uncertainties). The 
mean and median contributions of these pseudobulges to the 
total host galaxy light are 3.6% and 3.3%, respectively. 

The sole exception (J094112. 93+610340.7) has a 
marginal n = 2.0 ± 0.7, meaning we cannot say whether 
it is a classical bulge or a pseudobulge. Also, it is worth 
noting that the host galaxies with no detected pseudob- 
ulge are the highest-redshift sources in the sample; higher- 
resolution imaging could clarify both the status of this one 
exception and the pseudo/bulgeless nature of those sources 
with z > 0.06. 

It is therefore highly likely that most of these galaxies 
are truly bulgeless. However, we consider that all of the light 
from these components may be bulge light and consider it a 
robust upper limit on the contribution of a classical bulge. 
When no bulge is detected, we assume the upper limit to be 



tions (| Simmons fc Urrvll200Sl ). This additional uncertainty 
particularly affects the faint, compact second host galaxy 
components, but where a compact host component is de- 
tected we can nevertheless distinguish between bulge and 
pseudobulge in all but one system (described in Section 
[3]l ■ However, when calculating the upper limit to a possible 
bulge contribution to the galaxy (such as in Figure O, we in- 
clude all the light from even those compact host components 
firmly detected as pseudobulges. We therefore consider our 
bulge limits conservative upper limits. 



By construction, the sample is unambiguously disk- 
dominated. However, one source (J162511. 78+504202.1) ap- 
pears to be a merger of a galaxy with a strong bulge (n = 
2.65 + 0.18) and a companion, with tidal tails that resemble 
spiral arms. Another (J085903. 96+020503.9) is visually sim- 
ilar, but fitting indicates a disk-dominated (n = 0.6 ± 0.26) 
central component with a bulge-dominated companion. This 
system contains a broad-line AGN, and fitting the extended 
arms requires strong Fourier (asymmetric) modes. As the 
photometric redshift of the companion is consistent with a 
physical interaction between it and the primary source, this 
is likely a merger or post-merger and the arms may in fact 
be tidally induced features. Both galaxies are removed from 
our sample. 
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Figure 3. SDSS inverted r-band images and residuals from parametric fitting of AGN+hosts in Figure [2] Images are ordered by 
ascending redshift (as in Figure [2j; residual after subtracting best-fit model is below each image. Best-fit parameters are given in Table 

m 
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Table 1. Source properties for the 15 sources initially selected from Galaxy Zoo as potentially bulgeless AGN+host galaxies. After 
parametric morphological fitting, we retain the first 13 sources in the sample. All sources have an extended host galaxy component 
and nuclear unresolved point source. 11 of 13 sources also have a resolved, but compact, nuclear host component. Ten compact host 
components have n < 2 and are thus considered pseudobulgcs. The remaining nuclear host component has n ~ 2; within the uncertainties, 
we cannot distinguish whether this source has a classical bulge or a pseudobulge. The first source, J154832. 59+083516. 5, has two extended 
disks and a small pseudobulge. Most of the black hole masses are lower limits calculated by assuming the black hole is radiating at the 
Eddington limit; the three sources in the table with broad He? emission lines have firm black hole masses, but J085903. 96+020503. 9 is 
rejected from the sample due to a merging companion. 



2% of the host galaxy light. The minimum detected pseu- 
dobulge contribution is 1%; the maximum is 7.5%. 



3.2 Host Galaxy Properties 

The sample is insufficiently large to draw significant con- 
clusions about host galaxy properties, but it is worth not- 
ing that the galaxies, as seen in Figure are heterogenous 
examples of disk galaxies, comprising both barred and un- 
barred, red and blue, and tightly and loosely wound spi- 
ral samples. Even with this small sample, it is obvious that 
bulgeless galaxies hosting AGN are not restricted to a sin- 
gle morphological class. They span typical stellar mass and 
colour values for inactive late-type galaxies (Figure 2} . Al- 
though this sample does not contain host galaxies firmly 
on the red sequence, comparison to the inactive population 
is difficult owing to our emphasis on a pure rather than 
a complete sample. Disk galaxies on the red sequence are 
more likely to host growing bl ack holes than those in the 
blue cloud (|Masters et al.ll2010t ); further work is required to 
determine the fraction of bulgeless disk galaxies (both host- 
ing AGN and not) on the red sequence, as well as the overall 
population of bulgeless galaxies. 

Table [T] summarizes the results of parametric fitting. 
The Lhost/^AGN ratios range from 0.003 to 0.069. Extended 
Sersic profiles are disk-like (median n = 0.73), and the disk 
effective radii range from 2.4 to 15.2 kpc. 

One galaxy (J154832.59+083516.5) is best fit with two 
extended disks, which have similar axis ratios (b/a ~ 0.8) 
but are rotated by approximately 16° with respect to one 
another. The two disks are visually evident in the colour 
image (top left panel of Figure [2]| , and the fit improves with 
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Figure 4. u — r colour versus stellar mass for SDSS inac- 
tive late types with 0. 02 < z < 0.05 (black and grey contours; 
ISchawinski et al"1l2010h . with the 13 bulgeless AGN host galax- 
ies from this work shown in green. The bulgeless host galaxies 
span a range of masses and colours but are mainly located in the 
blue cloud at ma sses typical of the inactive disk galaxies from 
ISchawinski et all However, as the AGN host sample presented 
here is incomplete and was not selected in the same way as the 
inactive late-type sample, we caution against detailed compar- 
isons between these populations. 



the addition of the second extended disk component: the 
reduced goodness-of-fit parameter is xt = 1-032 with two 
disks versus 1.244 with a single disk (with ~ 40, 000 degrees 
of freedom), a significant improvement in the fit. 
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We note that the faint apparent companions to two 
of the galaxies shown in Figure [2] (J160534.64+323940.8 
and J110308. 04+072744. 4) are in fact projections of more 
distant background galaxies, according to photometric red- 
shifts. J160534. 64+323940. 8 also has an asymmetric spiral 
arm pattern, but asymmetric spiral features are not neces- 
sarily signatures of interaction, as disk instabilities can lead 
to bars, warps an d other asymmetric features wit hout ex- 
tern al interaction llSaha et al.ll2007l : ISellwoodll2010h . 

iBaldrv et al.l (|2006l ) derive a useful environmental mea- 
sure for SDSS galaxies closer than a redshift of 0.085. The 
local density for a galaxy is given by Eat = N/ {^:d 2 N ^j where 
djv is the projected distance to the Nth nearest neighbour 
that is more luminous than M r = 20, and E is determined by 
averaging the density determined using spectroscopic neigh- 
bours with that from using both photometric and spectro- 
scopic nei ghbours. We use the e xtension to SDSS DR6 de- 
scribed in Bamford et~aH l|2009h . 12 of our sample are in- 
cluded in the catalogue; values range from E = 0.076 to 
2.48 with a mean of 0.450, corresponding to the density of a 
typical field environment. For comparison, the entire SDSS 
has typical local density measures from 0.2 to 25. The envi- 
ronmental dependence of our systems is therefore markedl y 
different from that of merging galaxies (jDarg et al.ll2010h . 
supporting the idea that these bulgeless systems are free 

from recent merger activity. 

Stellar masses were derived following IBaldrv et al.l 

l|2006l ) using a best fit stellar mass-to-light ratio corrected 
for the observed dependence o nu-r colour. This ap proach, 
which was first introduced bv lBell fc de Jond (|200lh . is not 
as accurate as one based on full spectral fitting, but retains a 
simple relation between observed and derived quantities. It 
is also less likely to be distorted by the presence or absence 
of an AGN with luminosities like those in our sample. 



4 BLACK HOLE MASSES 

Two of the AGN in the sample (J120257.81+045045.0 and 
J133739. 87+390916.4) have broad emission lines. We used 
the Ha line width and flux to measure the black hole 
masses of these so urces, following iGreene fc Hoi l|2007h and 
Ijiang et al.l (|201ll ). We began by fitting the stellar contin- 
uum present in e ach spectrum using the GANDALF software 
ijSarzi et al.ll2006t >. Although GANDALF is designed to fit both 
absorption and emission lines, the profiles of the narrow 
emission lines in both objects are complex and asymmetric, 
so we instead subtracted the continuum fits from the data 
and analyzed the residual emission-line spectra ourselves. 
The process involved (a) modeling the profile of the [Oiii] 
5007 line as the sum of 2-3 Gaussian components, (b) fit- 
ting the [Oiii] model to the [Sn] 6716,31, [Nil] 6548,83, and 
narrow Ha lines, and (c) modeling the residual broad Ha 
components as a single Gaussian (J133739. 87+390916. 4) or 
Lorentzian (J120257.81+045045.0) profile. The results are 
displayed in Figure [5] For J133739. 87+390916.4, the broad 
Ha line has a width of 4950 km/s FWHM and a luminos- 
ity of 1.38 x 10 40 erg/s, which correspon d to a black-hole 
mass of 1.2 x 10 7 M Q (| Jiang et al.l 1201 ll ). The broad Ha 
line in J120257.81+045045.0 has a width and luminosity of 
2810 km/s and 1.9 x 10 40 erg/s, respectively, suggesting a 
black-hole mass of 4.2 x 10 6 Mq . A Gaussian fit to the broad 
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Figure 5. SDSS spectra (dark gray) of objects (a) 
J120257.81+045045.0 and (b) J133739.87+ 390916. 4, with spec- 
tral fits to narrow and broad components (blue). Narrow [Nil], 
[Sn], and Ha are fit separately in order to measure the isolated 
flux and width of the broad Ha component. 



Ha component of this object would imply a slightly larger 
black-hole mass, although it is inappropriate based on the 
spectral fits (Figure [5}. 

Lower limits on black hole masses for the remaining 
sample can be obtained using the bolometric luminosity of 
the sources. For sources like these where the host galaxy 
dominates the emission at most wavelengths, bolometric lu- 
minosities are typically obtained using corrections to either 
the hard X-ray or mid-infrared bands, where the AGN emis- 
sion dominates. All the sources are det ected by the Wide - 
field Infrared Survey Explorer ( WISE: IWright et alj|20ich 
with S/N > 5 . We use the waveleng th-dependent bolometric 
corrections of I Richards et all (|2006 ). which are not strongly 
dependent on AGN luminosity, to estimate the bolometric 
luminosity based on the longest-wavelength infrared data 
available for each source. Twelve of 13 AGN with bulge- 
less hosts are detected in the W4 band centered at 22 /im 
(Lbol ~ 10 x Z/22/jm); all are detected in the W3 band 
centered at 12 (Lho\ ~ 8 x Li2 f im). We verified that 
black hole mass limits calculated in this way are not sig- 
nificantly different from mass limits for this sample calcu- 
lated using a bolomet ric correction to the [Oiii] luminosity 
i|Heckman et alj|2004h . 

Given the bolometric luminosity, making the assump- 
tion that accretion is at the Eddington limit yields a lower 
limit on the black hole mass. Super-Eddington accretion 
is required for a black hole mass limit derived in this 
way to be overestimated; this is only rarely observed and, 
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when present, typica lly results in only a small change in 
observed luminosity (|CoIlin fc KawaguchJ|2004h . We note 
that the two sources with measured black hole masses 
are accreting at rates well below their Eddington limits 
(L bol /L E dd = 0.05 and 0.08 for J120257.81+045045.0 and 
J133739.87+390916.4, respectively). If these growth rates 
are typical of the sample, it implies the actual black hole 
masses of the sample are higher by at least ~ 1 dex than the 
computed lower limits. The lower limits obtained on black 
hole mass are substantial, ranging from ~ 10 5 to 10 6 M©. 
The black hole masses/mass limits are given in Table [1] 



5 DISCUSSION 

Massive pure disk galaxies hosting growing supermassive 
black holes offer a powerful probe of the co-evolution of 
galaxies and black holes in the absence of mergers (some- 
times called secular evolution). Their apparent rarity poses 
a significant challenge to using them to study secular means 
of black hole growth and arises from a combination of fac- 
tors. Firstly, in a universe where galaxy evolution proceeds 
via mergers, massive bulgeless disks formed without merg- 
ers a re expected to be very rare (e.g.. ISteinmetz fc Navarro] 
l2002h ; additionally, only a small fraction of those are ex- 
pected to host actively growing supermassive black holes 
at any given time. Identification of AGN via optical emis- 
sion lines is possible over volumes large enough to enclose 
a rare host galaxy sample, but such diagnostics miss opti- 
cally obscured AGN; the nuclear emission from those AGN 
that are detected can lead to visual confusion with a small 
bulge, impeding the identification of bulgeless galaxies host- 
ing optically-selected AGN. 

However, using a combination of Galaxy Zoo visual clas- 
sifications within the Sloan Digital Sky Survey and paramet- 
ric morphological analysis, we have selected 13 such galaxies, 
the largest sample of its kind to date. The unambiguously 
disk-dominated sample presented here was selected to favour 
purity rather than completeness and can thus be considered 
a robust initial sample of massive (M* > 10 10 M©) AGN 
host galaxies without classical bulges and with very small 
(or undetected) pseudobulges. 

The use of the Sersic index n = 2 to distinguish between 
dynamically cold pseudobulges, created by re-distributing 
disk matter in the absence of a merger, and dynamically 
hot (classical) bul ges, created in merge rs, is a dichotomy 
with little overlap jFisher fc Drory||2008l '). While the typical 
Sersic index f or a classical bulge is a function of luminos- 
ity/m ass fe.g-. lGraham fc Guzmanl2003l ; lGraham fc Worlevl 
120081 ) . classical bulges with n < 2 are very rare, especially for 
galaxies with comparable masses to this sample. It is thus 
highly unlikely that the compact host galaxy components 
with n < 2 are classical bulges. 

Additionally, the contribution of pseudobulges to these 
systems, derived from the fitting described above, is ex- 
tremely small, even compared to other samples of AGN 
host galaxies wit h histories thought to be dominated by 
secul ar evolution (|Orban de Xivrv et al . 2011; M athur et alj 
l201lh or to samples of in active massive bulgeless galax- 
ies |Kormendv et alj 20101) . Of the bulgeless giant galax- 
ies discussed in Kormendv et all only the extremely disk- 
dominated (3 of 19 in that sample) have comparable 



pseudobulge-to-disk ratios. In our sample, 12 of 13 host 
galaxies have either no second component or a small pseu- 
dobulge. Only one host galaxy has a compact component 
that is marginally consistent with a classical bulge, but this 
component is equally consistent with being a pseudobulge, 
contributing less than 4% of the total host galaxy light. 
The galaxies in the sample are thus fully consistent with 
being pure-disk AGN host galaxies, which places strong con- 
straints on their merger histories. 

In particular, these galaxies have no major mergers in 
their collective history, and their minor merger history is 
strongly limited. Simulations typically show that mergers 
with mass ratios greater than 1 : 10 for m a classical bulge 
jWalker et al.lll996l ; lHopkins et al.ll201ll ); the lack of bulges 
in this sample constrains its merger history to events with 
a smaller ratio (in the sense that the satellite is less than 
l/10th the main galaxy's mass). 

Note that some recent work suggests bulges may be sup- 
press ed in mergers up to a mass ratio of ~ 1 : 4 (I Brook et al.l 
|2012| ). which, if true, would allow a galaxy to remain bul- 
geless after merging with somewhat larger galaxies than 
the typical limit from most simulations. However, such a 
significant minor merger would leave signs of tidal debris; 
such signs could be detectable in images at the depths 
shown here, depending on the stage of relaxation (e.g., the 
m ore minor mergers in the spheroidal post-merger sample 
of ICarpineti et alJl2012T ). None of the galaxies in the sam- 
ple show evidence of tidal features at the SDSS depth, dis- 
favouring the notion that a minor merger event precipitated 
the obser ved black hole growth (deepe r imaging cou ld con- 
firm this; lKaviraill2010f ). Additionally, iBrook et al.l predict 
long-lived bar features in bulgeless galaxies that have under- 
gone a minor merger, but no excess of bars is observed in 
the sample, suggesting this is not a significant effect. Neither 
major nor minor mergers are driving the observed black hole 
growth. 

Given the bolometric luminosities of the sample, and 
assuming the bolometric luminosity Lboi (energy radiated) 
is related to the accretion rate rh (energy captured) by a 
radiative efficiency factor e, one can estimate the amount 
of matter falli ng onto each bla ck hole. Adopting a value 
of e = 0.15 jElvis et al.l l2002h yields accretion rates of 
0.003 < m < 0.03 M yr" 1 for the sample. This level is 
easily achievable with fe eding via dynamical ly cold accre- 
tion of minor satellites (|Crockett et al.ll201ll ). but it may 
also b e possible wi t h me rger-free processes alone. For ex- 
ample, ICiotti et alj (|l99ll ) calculated that stellar mass loss 
from a passively evolving stellar population (albeit in an el- 
lipti cal galaxy) could send m aterial toward a central SMBH 
Csee ICiotti fc Ostrikeri I2012J for a recent treatment, and 
iDavies et al.l 120071 for examples in local Seyfert galaxies). 
The amount of material that could be driven toward the 
nucleus within a disk galaxy remains uncertain and the 
details of such a process are unclear, but gas-rich disks 
with higher star formation rates and the increased potential 
for t ransfer of angular momentum via, e.g., bar instabili- 
ties (|Masters et alJl^OllT ) and viscous torques, present many 
purely secular opportunities for feedin g central sup ermassive 
black holes (for a detailed review, see IJogedbOOd l. 

Using black hole accretion rate and mass, one can esti- 
mate the time required to grow a seed SMBH to the observed 
mass. We assume the black holes have grown at the rate (i.e., 
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Figure 6. Black hole mass versus up p er lim its on bulge mass (left) and total host galaxy stellar mass (right). Host galaxy stellar masses 
are calculated following iBaldrv et al,l d2006h . Most of the sources have no detected classical bulge; we calculate robust upper limits to 
bulge masses (left panel; leftward arrows) using the pseudobulgc fraction of the r-band galaxy light profile. The single source with a 
disk-subtracted central galaxy profile that may be consistent with a classical bulge has a bulge mass indicated with a filled circle. AGN 
with broad-line emission have black hole masses calculated from Ha line width and flux (black squares). Vertical arrows indicate lower 
black hole mass limits calculated assuming radiation at the Eddington limit; supcr-Eddington accretion is req uired for a black hole to 
have a mass below the limit. The local relation between bulge/elliptical stellar mass and blac k hole mass fr om lHaring &: Rixl | |2004| ) is 
plotted as a dashed line; the shaded region indicates the observed scatter of 0.3 dex quoted bv lHaring fc Rixl , The black hole masses are 
higher than predicted by that relation based on the maximum bulge masses of these systems. However, all black hole masses/limits are 
consistent with the relation if the relation describes total stellar mass. If the Eddington ratios for the two systems with measured black 
hole masses are typical of the sample, the limits underestimate the black hole masses by at least ~ 1 dex, amplifying the offset in the 
left panel, but bringing the sample into good agreement with the correlation between black hole mass and total stellar mass. 



the luminosity) currently observed since the mass at which 
that rate was the Eddington rate, and assume Eddington- 
limited accretion for mass growth that occurred prior to that 
point. This effectively means the observed accretion rate is 
assumed to be the maximum rate (in M@ yr -1 ) at which 
the black hole has grown over its lifetime. Given the low 
rates observed, this is quite conservative, but it provides an 
estimate of the total time a black hole of the given mass 
would need to spend in an actively growing phase over its 
lifetime if the observed rate is typical. For the two systems 
with firm black hole masses, the time required is between 
~ 1 — 2 Gyr, depending on the seed mass (th e range given as- 
sumes seed masses be tween 10 2 and 10 5 ; I Volonteri et al.1 
l2008l : IVolonterill201Ch . 

Thus even if the currently observed accretion rate is 
the maximum rate over the lifetime of these AGN, the 
time required to grow SMBHs to the masses observed 
here is considerably less than a Hubble time. The SMBHs 
in bulgeless galaxies need only spend ~ 10% of their 
lifetimes in an actively growing phase, a similar fraction 
to that predicte d by independen t models and observed 
by others (e.g., Kauffman n et al. 20031; Hao et al. 2005; 
iHopkins fc HernquistfeOOrj ; iFiore et alj|2012h . This does not 
rule out the possibility of higher accretion rates in the past 



(e.g., with cold accretion flows at high redshift) and less 
time spent in an active growth phase, but such a phase is 
not necessary. Further characterisation of this requires firm 
black hole masses for the remainder of the sample. 

Whatever process has fed these SMBHs, it is clear that 
the formation of these very small pseudobulges does not cor- 
relate their properties with Mbh in the same way as that 
of classical bulges. Figure [5] compares black hole mass with 
the contribution to the stellar mass of the host galaxy from 
a bulge or pseudobulge component (under the assumption 
that the mass-to-light ratio of the compact galaxy compo- 
nent is equal to that of the whole galaxy) , for comparison to 
meas ured galaxy-black hole correlations. Such co rrelations 
(e.g. iMarconi fc Huntj 120031 ; lHaring fc Rixl |2004| ) compare 
black hole masses to classical bulge properties. However, 
only one object has a compact host galaxy component whose 
Sersic index is marginally consistent with that of a classical 
bulge. Because the rest of the galaxies in the sample have 
no detected classical bulge, for these sources the maximum 
contribution is calculated by using the mass of the pseu- 
dobulge as a strong upper limit on the mass of a classical 
bulge. Comparis on with the measure d correlation between 
these quantities (|Haxing fc Rixl [20041 ) shows that the black 
hole masses allowed by the derived limits are systemati- 
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cally above the measured correlation with bulge stellar mass. 
This result differs from som e previous studies of pseudob- 
ulges (e.g.. I Jiang et al 1 l201ll : iMathur et alll201lf l. however 
we note that the pure disk AGN host galaxies presented here 
are typically larger than those in previous studies, and/or 
the pseudobulges smaller. 

Figure [6] also shows the relationship between black 
hole mass and total stellar mass, for comparison to the 
same black hole-galaxy mass relation. Note that the rela- 
tion is based on elliptical/bulge masses for primarily bulge- 
dominated systems (26 of 30 galaxies in iHaxing fc Rixll2004l 
are classified as E or SO) , so Mbui ge ~ M«.tot for most of the 
systems on which the relation of Haring fe Rbj is based, in 
contrast to the total stellar masses of the bulgeless galax- 
ies in this sample. Because the disks in this sample have 
very different dynamical histories than elliptical or bulge- 
dominated galaxies, we do not expect their masses to cor- 
relate with black hole mass in the same way as bulges if 
different dynamical histories lead to different rates of black 
hole growth. 

However, the comparison of bulgeless galaxies to the 
correlation based on bulge-dominated galaxies shows very 
good agreement. One of the two systems in the sample for 
which absolute measures of black hole mass are available is 
consistent with the relation, while the other has a black hole 
mass just below the observed scatter (0.3 dex; but given the 
uncertainties in the stellar masses, it is marginally consis- 
tent). The systems with limits on the black hole masses are 
consistent with the relation, but given the Eddington ra- 
tios of the two systems with measured black hole masses, 
these limits are likely underestimates by at least 1 dex. If 
we assume the Eddington ratios of the two systems with 
measured masses are typical of the sample and apply them 
to the black hole mass estimates of the rem aining 11 sys- 
tems, 9 of the 1 1 fall within the scatter of the lHaring fc Rbd 
|2004 ) relation between black hole mass and total galaxy 
stellar mass (the remaining two are outside the scatter by 
approximately the same amount as J120257. 81+045045.0). 
That the black hole masses of pure disk galaxies may cor- 
relate with total stellar host galaxy mass in the same way 
as SMBHs in bulge-dominated and elliptical galaxies - de- 
spite very different galactic formation histories - indicates 
that the evolutionary processes driving the dynamical and 
morphological configuration of the galaxy stellar mass may 
not be fundamental to the growth of the central black hole. 

The results presented here should be read in the con- 
text of results from a growing number of simulations which 
show that the black hole-galaxy connection is a reflection 
of mutual correlations between these two components and 
the overall gravitational potential of the da rk matter halo 
l|Booth fc Schavdl20ld : IVolonteri et alj|201ll ). such that the 
black hole-galaxy rela tion is a natural outcom e of hierarchi- 
cal galaxy evolution jjahnke fc MacciollioTll ) regardless of 
the merger history of the galaxy. The observational evidence 
is less clear: some work ind icates a correlation b etween halo 
mass and black hole mass (|Bandara et al.ll2009l) . with some 
evidence that outliers to the M — a relation are not outliers 
on a similar relation between SMBH and dark matter halo 
|Bogdan et al.|[2~012T l , but others find no correlation between 
black hole mass and dark matter halo (|Kormendv fc Bender! 
l201lT) or between pseu dobulge or disk mass and BH mass 
(|Kormendv et al.ll20lj ). 



Regardless of whether the galactic-scale evolution of 
baryons is fundamental to the evolution of supermassive 
black holes, bulgeless galaxies hosting AGN provide a 
means of studying secular (i.e., merger- free/dynamically 
cold) pathways to supermassive black hole growth in rel- 
ative isolation compared to the majority of galaxies with a 
more complicated history of mergers and secular processes. 
They may eventually provide leverage to constrain the max- 
imum black hole growth possible via merger-free processes 
in all galaxies, and to separate the extent to which SMBH 
and galaxy properties correlate as a result of different evolu- 
tionary processes. The results presented here indicate that 
significant black hole growth (to at least ~ 10 7 Mq ) is possi- 
ble via pathways free of mergers, although a more complete 
treatment requires follow-up work to determine firm black 
hole masses for the remainder of the sample. 



6 CONCLUSIONS 

By using classifications of visual morphologies of Sloan Digi- 
tal Sky Survey galaxies, drawn from the Galaxy Zoo project, 
we have selected a large set of bulgeless face-on spiral galax- 
ies. A conservative initial selection identifies 13 of these 
galaxies which are unambiguously systems with growing 
black holes. Parameterized fitting of these galaxies provides 
stringent limits on bulge or pseudobulge mass, with the lat- 
ter typically contributing ~ 3% by mass on average. 

Two of the galaxies in the sample have broad-line AGN, 
and thus measurements of their black hole mass are possible. 
For the rest, infrared observations from the WISE mission 
allow us to place a lower limit on the black hole mass. The 
black hole masses are substantial, reaching ~ 10 7 Mq, and 
lie above those predicted by the local bulge-black hole mass 
relation, even when all the pseudobulge component is in- 
cluded as an upper limit to the mass of the classical bulge 
in each case. 

One the two black holes with measured masses is fully 
consistent with the relation between black hole and total 
stellar mass (the other is just outside the scatter). If the Ed- 
dington limits of the black holes with measured masses are 
typical of the full sample, 80% of the systems for which only 
lower limits are available have black hole masses consistent 
with predictions based on total galaxy stellar mass. Firm 
conclusions require further observations, but it is not incon- 
sistent with the idea that black hole mass is more closely 
related to the overall gravitational potential of the galaxy 
and its dark matter halo (which is dominated by the halo 
but traced by the total stellar mass) than to the dynamically 
hot bulge component. 

In any case, the presence of massive, growing super- 
massive black holes in bulgeless galaxies indicates that sec- 
ular evolution is an important part of the evolution of the 
galaxy population. Either significant black hole growth is 
possible even in the absence of significant mergers, or a dy- 
namical means to keep galaxies bulgeless despite mergers 
must be found. Future work will include an analysis of the 
more than 10,000 candidate bulgeless galaxies from which 
this sample was drawn in order to constrain the properties of 
this intriguing population; in particular, a search for bulge- 
less systems in mergers will distinguish between the two sce- 
narios left open. Observational follow-up of the small sam- 
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pie identified here, particularly in order to constrain more 
tightly the bulge properties and black hole masses is also ur- 
gently necessary. Although extending the work to higher red- 
shift will be challenging, Galaxy Zoo classifications for large 
Hubble Space Telescope studies hold the promise of identi- 
fying a similar set of galaxies out to a redshift of approxi- 
mately one. Even at low redshift, however, the systems iden- 
tified here present a stringent test of simulations of galaxy 
formation. 
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